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We show that any finite quantum system S can be coupled to a dephasing environment in such 
a way that the internal mechanism responsible for relaxation of observables acting on S can be 
effectively canceled. By adjusting this coupling, the difference between the initial and the long¬ 
time expectation values of any observable on S can be tuned to arbitrarily small however nonzero 
value. This statement is exemplified and visualized by numerical studies of relaxation in a generic 
one-dimensional system of interacting fermions. 
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I. INTRODUCTION 

There is ever-growing research in the field of dynamics 
of many-body quantum systems [I|. One of the central 
and still unsolved problems is related to equilibration of 
closed systems, with a particular emphasis on the relax¬ 
ation of observables relevant for experimental verification 
of various theoretical models. Moreover, in large-enough 
generic quantum systems one expects that equilibration 
is equivalent to thermalization. It is commonly accepted 
that the long-time averages of local observables in generic 
systems coincide with expectation values for the statisti¬ 
cal Gibbs ensemble M- However, our understanding of 
this process and the condition for its occurrence are far 
from being complete. For example, it is known that there 
is a relation between equilibration and integrability of the 
system Generally, nonintegrable systems are ex¬ 

pected to thermalize, while integrable systems do not ap¬ 
proach the Gibbs state but rather the generalized Gibbs 
ensemble (GGE) However, there are quantum 

systems which do thermalize despite being integrable pro¬ 
vided that they are prepared in certain states [l^- On 
the other hand, there are nonintegrable quantum systems 
which do not thermalize Q and the role of initial entan¬ 
glement between subsystems seems to be important. 

According to common intuition, relaxation of observ¬ 
ables is a hallmark of irreversibility [isj ]. If one consid¬ 
ers an expectation value of an arbitrary observable O 
and takes advantage of the spectral decomposition of the 
Hamiltonian into its eigenstates H\n) = En\n), one ob¬ 
tains for an initial state \ip) = C'ml™): 

{0)it)= C:Cra{n\0\m) 

+ ( 1 ) 

Invoking arguments adopted in the context of the eigen¬ 
state thermalization hypothesis Q , one expects that after 
sufficiently long time, due to destructive interference of 
oscillating terms in Eq. o, the observable can relax to 


its steady state value 

lim (0)(t) = lim — [ dt{0){t) 
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= Y. C:CMO\m). ( 2 ) 

Although finite closed quantum systems are strictly 
speaking periodic or quasi-periodic, the irreversibility is 
still manifested as low (or vanishing) probability of a pro¬ 
cess which reverses the relaxation. Moreover, if one at¬ 
taches the quantum system to an infinite environment, 
transforming a closed system into an open system, one 
expects that the relaxation becomes faster and ’more irre¬ 
versible’ due to the information loss or dissipation. This 
seems at least intuitively indisputable. However, in this 
paper we present a counter-example to this generally ac¬ 
cepted mechanism. We demonstrate that there exist en¬ 
vironments which cancel internal relaxation mechanisms 
of an arbitrary quantum system. More precisely, we show 
that any finite quantum system can be linearly coupled 
to a continuous set of harmonic oscillators in such a way 
that evolution of an arbitrary observable O satisfies: 

I lim {0){t) - lim(0)(t)| < e, (3) 

i— >oo t— s-O 

for arbitrary e > 0. In other words, the relaxation can¬ 
not be completely eliminated but can be made arbitrar¬ 
ily inefficient. This prediction is in direct contrast to Eq. 
([2|). The main idea behind this result is that a specially 
tailored environments may cancel the (destructively in¬ 
terfering) oscillations of the off diagonal matrix elements 
in Eq. O- Such environments belong to the class of 
purely dephasing environments [l8l| . They preserve the 
internal conservation laws of the quantum system, hence 
they are not generic or even typical. Applicability of 
the pure dephasing model has been discussed in vari¬ 
ous areas of quantum and mathematical physics 
Surprisingly, there are also real systems which can be ef- 
fectivel y d escribed by such a simple, highly symmetric, 
model [^. However, let us stress that pure decoher¬ 
ence remains credible only at the time scales significantly 


2 


smaller than the time scale relevant for exchanging en¬ 
ergy between system and its environment Hence, 

applicability of the pure dephasing under experimentally 
accessible conditions is usually at least disputable. 

The paper is organized as follows: In the next section 
we formulate and prove (in subsection A) the central re¬ 
sult of our work concerning the reversal of relaxation in 
dephasing environments. Further, in subsection B of Sec¬ 
tion II, we provide a simple example to illustrate reversal 
of relaxation in a generic system of interacting fermions. 
In Section III we present certain generalizations of the 
main result. Finally, we conclude our work in last sec¬ 
tion of the paper. 


II. MODEL AND MAIN RESULT 


We consider a quantum system, S, described by the 
Hamiltonian 


Hs = '^En\n){n\, (4) 

n 


with arbitrary (possibly degenerate) energy spectrum. 
The system is coupled to an environment, B, of non¬ 
interacting bosons 


Hb 


duj ioa^ {u!)a{oj), 


( 5 ) 


where the fields a{u}) satisfy [a(a;), a'l'(w')] = 6{uj — 
uj') [i^l- The details of the purely dephasing S-B cou¬ 
pling are given b y H s and Vb such that the total Hamil¬ 
tonian reads [la, 


H = Hs ® Ib +Ts ® Hb + Vs ®Vb, 


where 


Vb= dujg{uj) [a'‘'(a;) + a(a;)] , 
^0 

with a real valued g(oj) and 


= ^ln\n){n\. 
n 


( 6 ) 

( 7 ) 


( 8 ) 


Note that Hs and Vs commute. This particular prop¬ 
erty is a hallmark of a pure dephasing [l^. We assume 
also that initially the system and its environment are 
prepared in a separable state 


p(o)=i: Pnm|m)(n| (g) |H)(H|, (9) 

n,m 

where |H) is the bosonic vacuum. Towards the end of this 
paper we discuss more general class of the initial states. 


A. Proposition 

We consider long-time expectation values of local op¬ 
erators defined for the quantum system, O = Os ®Tb. 
If the system-environment coupling satisfies the propor¬ 
tionality relation 

7 m ~ 7 n O'- Em ~ En ( 10 ) 

then the difference | limt_,.oo (O) {t) — (O) (0) | can be tuned 
to arbitrary small value by an appropriate choice of the 
coupling function g{u}). Most importantly, a single tun¬ 
ing of g{uj) holds for all local operators O. 

For the pure dephasing one easily finds the time- 
dependent expectation value 

( 0 )(t) = E Pnm{n\0\m) 

O^nm it), ( 11 ) 

n,m 

where 


a„m(l) = (HI exp[it{En -I- ynUe + Hb)] 

X exp[-it{Em + ImVB + F1b)]|H). (12) 

Such simple result occurs due to the block-diagonal 
structure of the Hamiltonian, where each block describes 
a set of shifted harmonic oscillators: Em. + jmVB + Hb- 
An explicit form of a„m(l) has been derived/used many 
times in different contexts ranging from mathematical 
physics [ 2 ^ up to quantum information M- An explicit 
form of the amplitude anm{t) reads 

anmit) = 
where 


E{t) = 

1 

(M 

8_ 

(14) 


Jo ^ 


A{t) = 

poo 2 / N 

/ duj- —^ (1 — cos(a;t)). 

(15) 

Jo ^ 


Utilizing the square-integrability of q(uj)/u! one hnds 
from the Lebesgue-Riemann lemma [27| in the long-time 
regime (t 00 ) that 

E{t) —>■ Et and A(t) —>■ A = const (16) 


where 

noo 

E= dujg^{uj)uj~^ (17) 

^0 

nOO 

A = / dujg‘^{uj)uj~‘^ (18) 

Jo 

In order to reduce the effects of the internal relaxation 
one should tune the system-environment coupling in such 
a way that the oscillatory part in Eq. (fT^ drops out 

Em-E^- E{jI - 7 ') = 0 (19) 

while the exponential term ( 7 ^ — 7 n)^A remains small. 
In order to show that such a particular choice is indeed 
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possible, we introduce a cut-off frequency [l^, Wc, for 
the bosons in the environment and redefine the coupling 
function 


9 {^) = ( 20 ) 

Then, A oc u!~^ while E oc After eliminating 

the destructive interference of the off-diagonal matrix- 
elements, the exponential damping may become arbitrar¬ 
ily small, of the order of 1 /ujc 

log[amnit oo)] oc - ^ ^ 21 ) 

UJc 

While our general scheme does not require any particular 
form of g(uj) the results are most transparent for the 
standard parametrization of the coupling function [l 8 l| : 

9 \^) = exp(-c^M), (22) 

where p, > 0 corresponds to the (mathematically [ 2 ^) 
well behaving super-0hmic environment. In the latter 
case E = 1, while A = 


B. Example 

In order to exemplify our result with the help of a sim¬ 
ple but generic case, we study a one-dimensional system 
of L sites and Lj^ spin-less fermions with periodic bound¬ 
ary conditions |28l43ll | (i.e. a ring) with a Hamiltonian 
given by: 

L 

Hs = -th ^ |^exp(i()))c]+;^Cj -I- H.c. 
t=i 

L L 

-\-Ui 'y ( Tij^iTij + U 2 E nj+2nj (23) 

where rij = CjCj, th is the hopping integral, whereas Ui 
and U 2 describe first- and second-nearest-neighbor in¬ 
teractions, respectively. We take th as the energy unit 
th = 1, whereas time is expressed in units of h/th- We 
take also Ui = 1.4 and U 2 = 1. For such parameters 
Hs describes a generic metal characterized by a feature¬ 
less response to electromagnetic field and a normal diffu¬ 
sive transport. Moreover, such system thermalizes even 
when being decoupled from its surrounding [^. Numeri¬ 
cal studies have been carried out for L < 18. The reason 
behind introducing U 2 is to stay away from the integrable 
case, which shows anomalous relaxation [ij, [ij, 0 , |32j 
and charge transport [ssl - l^ . 

The mechanism of the reversal of relaxation (RR) holds 
either for all observables or for none. In the following, 
we show numerical results for a particle current 

L 

Os ='^i exp(i(^)c]+iCj -b H.C., (24) 


since {Os) vanishes whenever S is in equilibrium, hence 
large or even non-zero values of this observable imply 
that the system is in a non-thermal state. We con¬ 
sider a typical super-Ohmic environment with /r = 1 [see 
Ea. (l22l) ] when 


E{t) = t 


1 -b ^2^2 ’ 


m 


1 

Uc 1+ 


(25) 


While our qualitative conclusions do not depend on the 
initial state, we assume that the system is initially in 
a pure state Pnm = C'nC'm [see Eq. ©], where Cn = 
const > 0 if Re{n| 05 |m) > 0 and On = 0 otherwise. 
Such a state gives large initial current and has non-zero 
projection on large number (approximately one half) of 
energy eigenstates. Consequently, it leads to very clear 
numerical results for the relaxation of the particle cur¬ 
rent. 

In Fig. [T] we present the relaxation of current flow¬ 
ing in an isolated ring, i.e. in a closed system decoupled 
from any environment. Numerical results show that the 
fermionic system under consideration is generic and large 
enough so it relaxes (close) to equilibrium due to internal 
scattering processes even in the absence of any (dephas¬ 
ing) environment. In the same figure we show relaxation 
in the presence of a typical dephasing described by the 
following choice 7 ^ = [c.f. Eq.®], where is a 

random variable with a uniform distribution in an in¬ 
terval [0,.fo]- This case is very different from the tuned 
coupling described by Eo. lfTO)) . We notice that such a 
typical dephasing additionally accelerates the relaxation. 
Hence the model under consideration reproduces the well 
expected and intuitive results. 

The situation dramatically changes in the case of the 
fine-tuned coupling between the ring and its bosonic en¬ 
vironment satisfying Eq. cni). According to the Proposi¬ 
tion we expect reversal of the relaxation, i.e. after suffi¬ 
ciently long evolution time the expectation value of the 
current should approach its initial value. However, the 



Figure 1. (color online). Time dependence of the particle cur¬ 
rent in closed generic system (^o = 0 ) and for generic (ran¬ 
dom) coupling 7 ^ = ^rn, where is random variable with 
flat distribution in [0, ^ 0 ]. Here uJc = 1 has been assumed. 
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realistic systems are never perfect and one can expect 
that the condition in Ea. dlOl) is satisfied at most approx¬ 
imately. Let us assume that the the most optimal achiev¬ 
able tuning is given by 

7 m = ^rn ~ (26) 

where ^m S [ 0 ,^o] is again a uniformly distributed ran¬ 
dom variable. It describes a degree of quenched or frozen 
disorder present in our system due its imperfect prepa¬ 
ration. The time evolution of current flowing in the ring 
for different values of is presented in Fig. [21 In par¬ 
ticular, for an ideal case ^m = 0 condition in Eg. lfTUl) is 
satisfied and the RR clearly occurs. For small but non¬ 
vanishing values of there is still a wide time-window 
with significant degree of RR as indicated in Fig. |2l This 
nonmonotonic behavior with a wide plateau represents a 
hallmark of RR that could possibly be observed in simple 
quantum systems provided their couplings to the environ¬ 
ments could be appropriately tuned. 

Finally, we discuss the role of the bosonic character¬ 
istic (cut-off) frequency uJc- When this quantity is too 
small the coupling to the environment is not beneficial 
any more. Although RR is still possible, the off-diagonal 
matrix elements a„m decay in time due to the exponen¬ 
tial terms in Eq. (O. Results shown in Fig. 0 for a 
perfectly tuned coupling 7 ^ = Em — Eq indicate that a 
clear RR effect is observed already when tUc is one or¬ 
der of magnitude smaller that the typical energy scale of 
the quantum system, th- If Wc is of the same order of 
magnitude or even larger then the exponential damping 
becomes hardly visible. 

III. GENERALIZATIONS 

The assumptions behind the Proposition are not very 
restrictive, nevertheless they limit applicability of our 



Figure 2. (color online). Time dependence of the particle 
current for a partially tuned coupling 7 ^ = Em — Eq + ^m, 
where is random variable with flat distribution in [ 0 , 5 o]. 
uic = 0.1 has been assumed. The curves from the top (green) 
to the bottom (blue) are for ^0 =0, 0.01, 0.02, 0.03, 0.1, and 
1 , respectively 


main result to a certain class of problems. Here we dis¬ 
cuss which assumptions can be relaxed without a sig¬ 
nificant modification of Eq.(I3]). First, one may con¬ 
sider a more general class of initial separable states, 
Eq.(I 3 ), where bosonic subsystem is in the coherent state 
|S,^) = H(C)|f2), where D is the displacement opera¬ 
tor where C(w) is square-integrable and the prod¬ 

uct ((uj)g(uj) satisfies Lebesgue-Riemann lemma. This 
condition is satisfied, e.g. for a non-negative and square- 
integrable C(uj). Of course, for ^(w) = 0 one arrives back 
to the Proposition. 

It is known that the initial system-environment entan¬ 
glement can affect thermalization processes @ . A simple 
example of an entangled pure initial (not normalized) 
state is given by |n) (g) jS^;) -|- |m) <S> |Sp^) where |E!^) and 
|S(^) are noncolinear. Such a simple initial preparation 
of the pure dephasing models still allows us to find the 
exact time dependence of the total system-l-environment 
composite [^. Again, the reversal of relaxation can oc¬ 
cur provided that both C(a;)g(a;) and x(^)ff(‘^) satisfy 
the Lebesgue-Riemann lemma. 


IV. CONCLUSIONS 

The effect of an environment attached to a small sys¬ 
tem may occasionally be counter-intuitive or even un¬ 
predictable. Let us mention only two notable examples: 
stochastic resonance [4l| (both classical and quantum) 
and environment-induced entanglement . In our work 
we present another example of a counter-intuitive effect, 
which is the reversal of relaxation due to a coupling to 
the environment. We show that relaxation of an arbitrary 
observable acting on a finite quantum system can effec¬ 
tively be reversed by attaching the system to an infinite 
super-Ohmic dephasing reservoir. Since this mechanism 
requires a fine-tuning of the coupling between the energy 
levels and the bosonic bath, we expect that it could pos¬ 
sibly be realized in simple quantum systems rather than 



Figure 3. (color online). Time dependence of the particle 
current for a perfectly tuned coupling 7 ^ = Em — Eq and 
various cJc- 
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in complex generic cases. For the latter systems our find¬ 
ing will most probably remain only a matter of principle, 
in particular since such a coupling represents a non-local 
interaction. However, even for an imperfect tuning of 
the interaction between the system and the environment, 
one still finds a partial reversal of relaxation. This mech¬ 
anism may be realized in open quantum systems when 
the time scale of energy exchange with the environment 
is large when compared to other time scales of the sys¬ 
tem. Such an approximation has been effectively applied 
to describe quantum optical systems coupled to a sin¬ 
gle electromagnetic mode [l^. The effect studied in our 
work shows up as a broad time-window in which the ex¬ 
pectation values of observables are close to the values in 
the initial state. Even though the proposed mechanism 
for the reversal of relaxation is at this stage rather the¬ 


oretical, it may have an important impact in the area of 
quantum computing where one of the greatest challenges 
is controlling or removing quantum decoherence among 
interacting quantum systems. 

ACKNOWLEDGMENTS 

We acknowledges stimulating discussions with S. A. 
Trugman. This work has been supported by Pol¬ 
ish National Science Center under the grant DEC- 
2013/09/B/ST3/01659 (M.M. and J.D.) and by the Pl- 
0044 project of ARRS (J.B.). The work by J.B. was in 
part performed at the Center for Integrated Nanotech¬ 
nologies, a US Department of Energy, Office of Basic En¬ 
ergy Sciences user facility. 


[1] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- [22] 
tore, Rev. Mod. Phys. 83, 863 (2011) 

[2] S. Goldstein, J. L. Lebowitz, R. Tumulka, and N. Zanghi, [23] 
Phys. Rev. Lett. 96, 050403 (2006) 

[3] M. Srednicki, Phys. Rev. E 50, 888 (1994) 

[4] N. Linden, S. Popescu, A. J. Short, and A. Winter, 

Phys. Rev. E 79, 061103 (2009), [24] 

[5] C. Gogolin, M. P. Muller, and J. Eisert, 

Phys. Rev. Lett. 106, 040401 (2011) 

[6] A. Riera, C. Gogolin, and J. Eisert, 

Phys. Rev. Lett. 108, 080402 (2012) [25] 

[7] S. R. Manmana, S. Wessel, R. M. Noack, and A. Mura- [26] 

matsu, Phys. Rev. Lett. 98, 210405 (2007) [27] 

[8] L. F. Santos, A. Polkovnikov, and M. Rigol, 

Phys. Rev. Lett. 107, 040601 (2011) [28] 

[9] M. Mierzejewski, T. Prosen, D. Grivelli, and 

P. Prelovsek, Phys. Rev. Lett. 110, 200602 (2013), [29] 

[10] J. Larson, Journal of Physics B: Atomic, Molecular and Optical 

[11] M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, [30] 
Phys. Rev. Lett. 98, 050405 (2007) 

[12] D. M. Kennes and V. Meden, [31] 

Phys. Rev. B 82, 085109 (2010) 

[13] M. Kollar, F. A. Wolf, and M. Eckstein, [32] 
Phys. Rev. B 84, 054304 (2011) 

[14] A. C. Cassidy, C. W. Clark, and M. Rigol, [33] 

Phys. Rev. Lett. 106, 140405 (2011) 

[15] M. Fagotti, M. Collura, F. H. L. Essler, and P. Calabrese, [34] 
Phys. Rev. B 89, 125101 (2014) 

[16] M. Mierzejewski, P. Prelovsek, and T. Prosen, [35] 

Phys. Rev. Lett. 113, 020602 (2014) [36] 

[17] M. Rigol and M. Srednicki, [37] 

Phys. Rev. Lett. 108, 110601 (2012) 

[18] H.-P. Breuer and F. Petruccione, The Theory od Open [38] 
Quantum Systems (Oxford University Press, 2003). 

[19] J. Dajka, M. Mierzejewski, and J. Luczka, Phys. Rev. A [39] 
77, 042316 (2008). 

[20] K. Roszak and P. Machnikowski, Phys. Rev. A 73, [40] 

022313 (2006). [41] 

[21] J. H. Reina, L. Quiroga, and N. F. Johnson, Phys. Rev. 

A 65, 032326 (2002). [42] 


J. Dajka and J. Luczka, Journal of Physics A: Mathe¬ 
matical and Theoretical 41, 442001 (2008). 

D. I. Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, 
J. M. Gambetta, A. Blais, L. Frunzio, J. Majer, B. John¬ 
son, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, 
Nature 445, 515 (2007). 

O. Bratteli and D. W. Robinson, Operator Algebras 
and Quantum Statistical Mechanics: Equilibrium States. 
Models in Quantum Statistical Mechanics (Springer, 
2003). 

J. Dajka and J. Luczka, Phys. Rev. A 77, 062303 (2008) 

R. Alicki, Open Syst. Inf. Dyn. 11, 53 (2004). 

S. Bochner and K. Chandrasekharan, Fourier transforms 
(Princeton University Press, 1949). 

M. Biittiker, Y. Imry, and R. Landauer, 

Physics Letters A 96, 365 (1983) 

G. Blatter and D. A. Browne, 
Iffiityksltte,^. S4S)7b3(®a:|>988) 

R. Hiibner and R. Graham, 

Phys. Rev. B 53, 4870 (1996) 

R. Landauer and M. Biittiker, 

Phys. Rev. Lett. 54, 2049 (1985) 

R. Steinigeweg, J. Herbrych, P. Prelovsek, and 
M. Mierzejewski, Phys. Rev. B 85, 214409 (2012) 

M. Mierzejewski and P. Prelovsek, 
Phys. Rev. Lett. 105, 186405 (2010) 

M. Mierzejewski, J. Bonca, and P. Prelovsek, 
Phys. Rev. Lett. 107, 126601 (2011) 

T. Prosen, Phys. Rev. Lett. 107, 137201 (2011) 

M. Znidaric, Phys. Rev. Lett. 106, 220601 (2011) 

J. Sirker, R. G. Pereira, and I. Affleck, 

Phys. Rev. Lett. 103, 216602 (2009) 

R. Steinigeweg and W. Brenig, 

Phys. Rev. Lett. 107, 250602 (2011) 

J. Dajka, M. Mierzejewski, and J. Luczka, 

Phys. Rev. A 79, 012104 (2009) 

J. Dajka and J. Luczka, Phys. Rev. A 82, 012341 (2010) 
L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, 
Rev. Mod. Phys. 70, 223 (1998) 

F. Benatti, R. Floreanini, and M. Piani, 
Phys. Rev. Lett. 91, 070402 (2003) 





